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We present an approach to the computation of the non-Abelian statistics of quasiholes in quantum 
Hall states, such as the Pfaffian state, whose wavefunctions are related to the conformal blocks of 
minimal model conformal field theories. We use the Coulomb gas construction of these conformal 
field theories to formulate a plasma analogy for the quantum Hall states. A number of properties of 
the Pfaffian state follow immediately, including the Berry phases which demonstrate the quasiholes' 
fractional charge, the abelian statistics of the two-quasihole state, and equal-time ground state 
correlation functions. The non-Abelian statistics of multi-quasihole states follows from an additional 
assumption. 

Introduction. Recent investigations of the possibility of finding non-Abelian braiding statistics in the quantum Hall 
regime have used a remarkable observation due to Moore and Read as a springboard They drew attention 

to the fact that certain wavefunctions which have been proposed to describe incompressible quantum Hall states 
could be represented as conformal blocks in corresponding conformal field theories. This is not a statement about 
dynamics in the quantum Hall regime^], but merely a statement that some wavefunctions could be reconstructed 
from conformal field theories. Hence, it might not appear particularly useful. However, this observation suggests 
a way around a technical impasse which impedes the study of non-Abelian quantum Hall states. The problem and 
a strategy for its circumvention noted by Moore and Read may be summarized as follows. Suppose we initially 
eschew any attempt to realistically account for the physics of the long-range Coulomb interaction, finite thickness 
of the quantum well, inter-Landau-level mixing, and other complications, and, instead, work with an unrealistic but 
tractable model Hamiltonian whose exact zero-energy eigenstates (ground state and multi-quasihole states) can be 
explicitly constructed. The statistics of these quasiholes may be determined from the exact eigenstates using the 
Berry phase technique || . If the model Hamiltonian may be adiabatically connected to the one which governs the 
actual physical 2D electron gas - i.e. if the two Hamiltonians are in the same universality class - the statistics of the 
quasiholes and other 'topological quantum numbers' will be the same for both. In such a case, the calculation of 
these properties of the idealized Hamiltonian will be relevant to the real system. The stumbling block is that Berry 
phases (or Berry matrices in the non-Abelian case) are determined from matrix elements of the eigenstates and further 
progress is impossible unless we can find a method for calculating the matrix elements of the correlated many-body 
wavefunctions of interest. 

In his pioneering work on the fractional quantum Hall effect ]j| , Laughlin took advantage of the equality between 
the squared modulus of his wavefunctions and the Boltzmann weight of a classical two-dimensional plasma. The 
conventional wisdom on plasmas implies that the Laughlin states are quantum liquids || and has also facilitated 
calculations of off-diagonal long-range-order in these states |jj . The original calculation of the fractional statistics of 
the Laughlin quasiholes |(| relied on the incompressible liquid nature of the states but did not explicitly utilize the 
plasma analogy. However, the calculation may be rephrased 0] in the following way so that its dependence on the 
plasma analogy is emphasized. Let us write the wavefunction for a Laughlin state with two quasiholes at r\\ and 772 
in the following form (£0 is the magnetic length): 

k i>j 

Ordinarily, we would choose the wavefunction to be a single-valued function of the quasihole coordinates (which are 
simply some parameters in an electron wavefunction) by, for instance, replacing the first factor by its modulus or 



At least on the face of it. In fact, in many cases the conformal field theory which reproduces the wavefunctions is the same 
as the one which describes the gapless excitations at the edge. If the effective field theory of the bulk state is a Chern-Simons 
theory, this relationship can be understood in fight of Witten's discoveries |lo) : the states of the Chern-Simons theory which 
describes the bulk are isomorphic to the conformal blocks of the CFT which reproduces the bulk wavefunction. At the same 
time, the boundary theory which is concomitant with the Chern-Simons theory is precisely the CFT of the edge. The equality 
of these CFTs follows from the general covariance of the Chern-Simons theory. 
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removing it altogether. In (Q), we have intentionally chosen the phase of the wavefunction such that it is a multi- 
valued function of the quasihole coordinates (but, of course, single- valued in the electron coordinates); the net effect 
of braiding the two quasiholes will be the Berry holonomy together with the explicit monodromy of ([!]) . Only this 
combination is physically observable, and by changing the phase of the wavefunction - i.e. performing a gauge 
transformation - we can divide this phase between the Berry holonomy and the explicit monodromy in any way we 
wish. As we will see momentarily, the choice ([!]) is particularly convenient because it is the 'fractional statistics gauge' 
|)| in which the Berry holonomy vanishes and the explicit monodromy is the whole story. A second advantage of (|l|) 
is that it is a holomorphic function of 771 , 772 up to the Gaussian factor. Hence, the Berry holonomy, 7, acquired when 
rji is taken adiabatically around 772 is: 



where $ is the magnetic flux enclosed within the path taken by 771 and $ is the flux quantum. The second equal sign 
follows from the fact that {^^^2 1 depends only on 771 and not on 771, except through the Gaussian factor, which yields 
the second term on the right hand side. This term represents the fractional charge of the quasiparticles and does not 
depend on whether or not 771 encircles 772 . We now argue that as a result of the prefactor chosen in (0) , the first term 
vanishes. This follows from the observation that this inner product is the partition function of a classical plasma of 
charge y/m particles at temperature T = 1 with a neutralizing background of density l/2w\/m and two fixed charge 
1/y/m particles at 77! and 772 : 

O^Wl^iW = / l[ d Z i e " ° ( 3 ) 

(Usually, one takes T — m and has charge m electrons and charge 1 quasiholes, but the above convention is more 
natural for what follows.) When 1 771 — 772 1 is larger than the Debye screening length, the interaction between these 
two charges will be screened, so the partition function will be independent of 771, 772. In other words, (\E , - r7lTJ:2 |^E r T 7l972 } 
is independent of 77! and 772 so long as they are well-separated. Therefore, the Berry holonomy vanishes and the 
fractional statistics is given by the explicit monodromy of ([l]). 

Unfortunately, this argument appears limited to the Laughlin states; the squared moduli of the states which are 
promising hunting grounds for non-Abelian statistics, such as the Pfaffian state, 

*Pf = Pf (— IT" ) Ufa-**) 3 e"^ E ' Zl ' 2 (4) 

V %i Zj J . 

(where Pf (^ z .^ z . J — A y Zl l Z2 Z3 - Z4 is t ne antisymmetrized product over pairs of electrons) do not seem to be 
equal to the Boltzmann weight of a classical plasma, at least not in the most obvious way. A further complication arises 
from the fact that there is a degenerate set of multi-quasihole states in the putative non-Abelian cases, so off-diagonal 
matrix elements must also be calculated in order to obtain a Berry matrix. It is at this seemingly hopeless point 
that Moore and Read's observation can come to the rescue. Recall that the conformal blocks of certain conformal 
field theories are related to many of the quantum Hall states - including the Laughlin states and the Pfaffian state 
- which are exact ground states of model Hamiltonians. It is natural to conjecture that the quasihole statistics are 
equal to the monodromy properties of the conformal blocks (see also Q). These statistics are encapsulated in the 
concomitant Chern-Simons theories which describe the long- wavelength physics of the quantum Hall state [Q. 
Some justification for this conjecture, on rather general grounds, was given in ||, where it was used to conclude that 
the braiding matrices of the 2?7-quasihole states in the Pfaffian are embedded in the spinor representation of SO(2n). 

In this paper, we give further arguments (which we believe can serve as an outline for a proof) supporting this 
conjecture. We formulate a plasma analogy for non-Abelian quantum Hall states which makes use of the fortunate fact 
that all of the minimal model conformal field theories have a Coulomb gas description. At a first glance, this seems 
tantalizing but insufficient for our purposes because the construction involves complicated contour integrals of the 
locations of screening charges. However, following Mathur |pj| , we rewrite these contour integrals as two-dimensional 
integrals. As a result, we find that sums of the desired inner products can be expressed as the partition functions of 
plasmas with an additional large number of screening charges (which only enhance the screening property which we 
invoke) . The individual inner products are obtained by analytic continuation. This approach is applied to the Pfaffian 
state, which is associated with the c = 5 + 1 conformal field theory. Our plasma analogy enables us to compute the 
electron density matrix, the charge of the quasihole excitations (which, we argue, have localized charge distributions), 
and the Abelian statistics of the two-quasihole state. With one assumption regarding the analytic continuation, wc 
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can show that the Berry matrices of states with an arbitrary number of quasiholes are trivial and the non-Abelian 
statistics can be obtained entirely from the monodromy properties of the conformal blocks, in agreement with Q. 

Background Charge Construction and a Plasma Analogy. Let us momentarily take one step back before taking two 
steps forward. Recall that the Laughlin state can be expressed as a conformal block of the c = 1 theory, which is the 
theory of a free bosonic field <p: 



1 

2^ 



dipdtp (5) 



The Laughlin state is given by the conformal block of the vertex operator V^^(z) —: exp(i^/mipfi(z)) :, (tpn is the 
holomorphic part of the boson field) which plays the role of the electron, together with the operator for a neutralizing 
background charge. 

*l/m = ( e iv/ ™^ Zl ) e iv/ ™*^ 2 ^ e Vm0(zjv) e ~i J d 2 z ^h~p 4>(z)^ ^ 

Conformal blocks with V 1 /^i(r] i ) =: exp{iipn{r}i)/^/m) : are equal to states with quasiholes at The product 
of the Laughlin wave function and its complex conjugate can be written as the product of a holomorphic and an 
antiholomorphic conformal block or, simply, as a correlation function of the left-right symmetric field V y /j^(z i z) =: 
exp(i v / m((/?fl(z) + ipi,{z))) :. These correlation functions admit the Coulomb gas or plasma description of (||). The 
Pfaffian state, however, involves in its description a c = ~ conformal field theory. Can we describe it, too, by means 
of a Coulomb gas? 

Naively, the answer is yes since any conformal field theory admits a Coulomb gas description, as was discovered by 
Feigin and Fuchs and Dotsenko and Fateev (see JL2| and references therein). The Dotsenko-Fateev picture involves 
putting a certain background charge at infinity and introducing screening operators in a c = 1 conformal block. Then 
any conformal block of any conformal field theory can be expressed as 

F(z 11 z 2 ,...)=j c n^^no,K)n^.(z fc )^ (7) 

where the vertex operators O(wj) are the screening operators and Vk(zk) are the vertex operators of interest. The 
integration of the screening operator locations is along certain contours, Ci. The choice of contours determines which 
conformal blocks we obtain. Notice that the block on the right hand side of (0) is a standard Coulomb gas. Therefore 
it's tempting to conclude that we have a plasma analogy for any wave function connected to conformal field theory. 

However, this conclusion is premature. The contours for the integrals of (Q) are in general hard to define. To compute 
inner products we would have to introduce holomorphic and antiholomorphic screening charges, each integrated over 
its own contour. Since, in general, the number of screening charges we need will be of the order of the number of 
the electrons, it is unclear how we can compute the desired inner products. It is also unclear if we are allowed to 
exchange the order of integration over the positions z% of the electrons (which we will eventually have to do in order 
to compute the inner products) with the contour integrations over the screening charge locations since these contours 
are supposed to wind around the electron positions in a complicated way p^ |. 

A way out of this quagmire can be found in a paper by Mathur O]. He suggested the following. If, instead of 
considering holomorphic vertex operators, we consider full vertex operators which are the products of the holomorphic 
and antiholomorphic ones, V(z, z) = V(z)V(z), then it can be shown that 

j J]dV(n°K'%)II F ( Zfc ' Sfc )/ =^C n F n {z 1 ,z 2 ,...)F n {zi,z i ,...) (8) 

•* i \ j k In 

where the sum over n is the sum over different conformal blocks F n (different choices of contours in (]?])) while the 
C n are the so-called structure coefficients, which make the combination on the right-hand-side of (|^) a single valued 
function. 

The advantage of (||) over (0) is that the integration of the screening charges is over the entire complex plane, 
thereby putting the screening charges and electrons on equal footing in the computation of inner products. The 
disadvantage of (Jfy is that instead of producing the squared modulus of a particular conformal block (corresponding 
to a particular wave function) , it gives us a linear combination of them. We will discuss ways of getting around this 
problem later. For now, let us concentrate on cases in which there is only one conformal block, making the summation 
over n in (^) unnecessary. 

As an example, we consider the Pfaffian state (||) , which is a conformal block of the c = | + 1 theory |l],|j . The 
Pfaffian factor is the c = \ part, while the rest comes from the c = 1 part and has the standard plasma form of a 
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Laughlin state. The Laughlin plasma consists of charge \[2 electrons. According to [O], the c = i part can also 
be described as a plasma with N charges of magnitude 3/\/6 or —2/y/& associated with the electrons together with 
rs Y screening charges of magnitude 4/\/6 or — 3/v6 which also roam the entire complex plan and the background 
charge — l/v6 which has to be taken to the infinitely far away point. In other words, the background charge, being 
at infinity, does not interact with other charges, but shifts the tolal charge balance. The sum of all the other charges 
of our plasma has to be equal to minus the background charge (details can be found in Jl2| ) . We note that at large 
N it does not affect the neutrality of the plasma. 

We note that each operator can be represented by either of two different charges; this is a characteristic feature 
of the background charge construction. We have a certain freedom in how we pick either of these charges for each 
operator, and there is more than one equivalent representation of each conformal block. One particular construction 
is discussed in |l3| . However, these details are not relevant at the moment. For our purposes, it is sufficient to know 
that, together with the background charge, our plasma is neutral, and that, roughly speaking, we need one screening 
charge for each singularity encountered in the wave function. There are y singularities in the c — \ part of (Q), 
making the total number of screening charg GS OC "2" ; cL number which is large when the total number of the electrons 
N is large. 

Therefore we arrive at the central point of this paper: the Pfafhan state can be mapped to a plasma with three 
species of charges which interact via two different logarithmic interactions, i.e. 

|* Pf | 2 = f Y[d 2 w A e -m*i,n a ,wA) (9) 

3 A 

The inverse temperature, /3, is 1, and $ is the potential energy of a classical two-dimensional gas of logarithmically- 
interacting particles corresponding to: 

(a) electrons, which have charge \[2 with respect to the first interaction (c = 1) and charge 3/V§ or —2/^/6 with 
respect to the second (c = 1/2) 

(6) quasiholes, which have charge 1/2^/2 with respect to the first interaction (c = 1) and charge 3/(2\/6) or — 1/2\/6 
with respect to the second (c = 1/2) 

(c) the screening charges, which are neutral with respect to the first interaction (c = 1) and carry charge 4/-\/6 or 
— 3/y6 with respect to the other (c = 1/2) 

(d) the background charge, which is also neutral with respect to the first interaction (c — 1) and carry charge 
— 1/V6 with respect to the second (c = 1/2). 

The squared norm of the Pfaffian state is the partition function of this plasma. The squared modulus of the Pfaffian 
wavefunction (^|) - which, naively, is a mess - is, in fact, the Boltzmann weight of this plasma with the screening 
charges already integrated out. The integration of the N/2 wa's does not affect the c = 1 interaction and, at least 
partially, screens the c = 1/2 interactions among the z^'s and rja's. Hence, if we write @ as the Boltzmann weight of 
some effective Hamiltonian, <I>, it will have two species of particles which interact logarithmically through the c = 1 
interaction and also through some weaker interaction (perhaps short-ranged) which is the partially screened c = 1/2 
interaction. Though this latter interaction may be very complicated, we expect that (^) still essentially describes a 
plasma even after the w^'s have been integrated out. 

Quantum Liquid Nature and Electron Density Matrix. This observation allows us to immediately compute various 
properties of this state. Following Laughlin, we anticipate that it is a liquid because the associated plasma is above 
the crystallization temperature. As a result, we can compute its electron density matrix 

N 

P(z, A= / II d 2 z^* Pf (z', z 2 , . . .)*p/(*, z 2 , . . .) (10) 

^ i=2 

We note that p{z, z) is given by a partition function of our plasma with one charge held fixed at z\ in a liquid, it must 
be independent of z. On the other hand, 

p(z, z') = exp (- i(|z| 2 + |z'| 2 )) G(z, z') (11) 

where G{z, z') is analytic in its arguments. The only way to satisfy these two constraints is if 

p(z,z') = exp(-i(|z| 2 + |z'| 2 ) + izz') (12) 
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just as in the Laughlin states. There is no sign of pairing or any delicate power-law correlations in the electron 
Green function. Indeed, we don't expect such correlations in multi-electron Green functions either; all of these should 
be Gaussian or exponentially localized as a result of the gap. The interesting correlations should appear in Green 
functions of the non-local order parameter 0]. These correlations are related to the 'topological ordering' of the state 
and, hence, should be more stable against local perturbations than putative subtle correlations in the electron Green 
functions would be. Unfortunately the plasma analogy in the form presented here does not seem to allow to compute 
them. 

Finally, we can argue, following Laughlin ||, that the charges of the fractionally charged particles are localized 
around their positions, due to the effective plasma repulsion. The quasiholes do not have power-law tails in their 
charge distribution as one might have naively guessed from the form of the wavefunctions and from the fact that they 
carry half of a flux quantum. 

Fractional Charge and Abelian Statistics in the Two-Quasihole State. There is a unique two-quasihole state (see Q]) 
if the quasiholes are fixed at positions 771 and 772. It equal to a conformal block with twist fields and appropriate c = 1 
vertex operators inserted to represent the quasiholes [QJi) . According to our plasma analogy, the squared norm of the 
two-quasihole state is the partition function of our two-component plasma with two external charges (of magnitudes 
l/2\/2 with respect to the c = 1 interaction and 3/(2v6) or — 1/2 v6 with respect to the c — 1/2 interaction) held 
fixed at 771 and 772- The plasma will screen these charges, so the partition function will be independent of 771 and 772 
when their separation is larger than the screening length: 

*p7 hOl0S (^1^2)|^p7 h ° 1CS (?7l, m) ) = Ztwocomp. pW^l,^) = COllSt. (13) 

and consequently the non-trivial part of the Berry connection of these wave functions vanishes (the trivial part, coming 
from the Gaussian, gives the fractional charge of the quasiholes, 1/4) 



iTr2 — holes/ m \ 



d 



*pj holes (77i, 772) ) — term coming from the Gaussian = (14) 



_d_ 

di]i 



* P 7 hol0S (r/i, ? 72)|*p7 holcs (r/ 1 , 772) ) = 



Therefore we can read the fractional statistics of the quasiholes of the Pfaffian state from the monodromy of the 
conformal block which, is this case, is determined by a factor of (771 — 772)°. Interestingly, a = or, in other words, 
in a two-quasihole state, the two quasiholes are bosonic with respect to each other [pM]. The fractional statistics 
which is concomitant with the fractional charge is cancelled by the fractional statistics resulting from the Pfaffian 
part of the wavefunction. Since the charge and statistics of the quasiholes are different, it should actually be easier 
to experimentally disentangle them in the Pfaffian state than in a Laughlin state. 

Non- Abelian Statistics of 2n-Quasihole States. Now let us take up the less trivial case of four quasiholes. There 
are two possible four-quasihole wave functions when the quasiholes are located at 771, 772, 773, 774. We call them and 
^1 , following 0. According to (0), the plasma analogy allows us to compute only the combination 



(* |*o) + ( / = Ztwocomp.pla S ma(?7l, f?2,7i3, 7/4) = COnst. (15) 

This is not enough to argue that the braiding statistics of the quasiholes can be read off the conformal blocks. What 
we need is the somewhat stronger statement, 

(* i |* i )=C ii , i,i = 0, i (16) 

where is some constant (independent of the positions of quasiholes) matrix, which, if true, implies that the 
nonabelian Berry connection vanishes 



A i3 = ( 



d_ 

drj 



*./) = |^=0 (17) 



To derive ( |l6| ) we apparently need to generalize (||). However, we can instead look at the problem in the following 
way. (||) represents a single valued function, understood as a function of 77^ and 77^ = 77*. However, as a function of r\i 
only, with fji viewed formally as an independent variable, (||) is not single- valued. Moreover, by taking one r\i around 
another (while keeping fji at bay) we can transform (||) into an arbitrary combination of conformal blocks. For example, 
if we take 771 around 772 while holding 771 and 772 as well as 773, 773, 774, 774 fixed, then ty — * *o while \&i — * — If 
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we instead take r\\ around 7/4, while holding the rest fixed, "Jo — > and \& 1 — > "Jo- Hence, with these two braiding 

operations we can transform (* |*o) + ( * 1 1* \ ) into (* |*o) - i|, and ^* |*i \ ± ^*j|^o/ 5 from 

which we can obtain the desired inner products (|l6|). In other words, after performing these braiding operations, 
we can express the above combinations of inner products as analytic continuations of ^two comp. plasma (jli > i]2iV3 5^74)- 
The crux of the matter is now: can these analytically continued partition functions also be interpreted as partition 
functions of plasmas, albeit with modified interactions? If so, then we can, as usual, invoke plasma screening and 
conclude that all of the inner products are independent of the positions of the quasiholes. 

At the moment, we do not have a proof that this is so. However, the following line of reasoning strongly supports this 
contention. Earlier, we argued that integration of the wa's should result in some sort of plasma with partially screened 
interactions. If we then proceed and integrate out the Zi's, we expect total screening of the interactions between the 
rjaS and, hence, the result (|l5|). Suppose, instead, that before integrating the Zi's we first carry out one of the analytic 
continuations described above. This will not effect the c = 1 interaction since the area enclosed by the path of the 
r] a 's can be made arbitrarily small and the c = 1 interaction will only be sensitive to the analytic continuation if one 
of the Zj's crosses this path. The c = 1/2 interaction can be affected, but only through the introduction of phases 
which can change as Zj's approach rj^s. If these phases interfere destructively, the off-diagonal matrix elements vanish 
and the diagonal matrix elements are equal and constant. At worst, they can interfere constructively, in which case 
the other matrix elements are in the same situation as (U). This scenario, too, leads to the desired result ([r]) if we 
again appeal to plasma screening. 

If this argument proves to be correct, we can conclude (by straightforwardly extending it to the 2n-quasihole case) 
that the statistics of quasiholes in the Pfaffian state is given by spinor representations of SO(2ri), as was conjectured 
in§. 

Discussion. As we discussed in the introduction, some quantum Hall wavefunctions can be reproduced by conformal 
field theories, and it is natural to conjecture that the braiding properties of the former are given by those of the latter 
In this paper, we have attempted to answer the question 'why?' The basic reason, we submit, is that these 
states are all related - through the background charge construction (l2) - to Coulomb gases. As a result of plasma 
screening in these Coulomb gases, wavefunctions given by conformal blocks have trivial Berry connections; their 
braiding properties are determined entirely by the monodromies of the conformal blocks. 

Thus far, we have focussed on the Pfaffian state. However, it is clear from our discussion that the arguments we have 
used are quite general and may be applied to any quantum Hall state with wavefunctions given by the conformal blocks 
of a theory with a background charge construction. To the casual observer, the background charge representation of 
a conformal block may seem like a complicated integral represention. However, it is perfectly tailored for the present 
context because the electron coordinates must also be integrated in the calculation of Berry matrices. Consequently, 
the screening charges are on the same footing as the electrons and the representation (JsJ) is a natural generalization 
of the Coulomb gas representation of the Laughlin states. 
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